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ABSTRACT 

We apply the stellar population synthesis code by Cid Fernandes et al. to model the 
stellar contribution for a sample of 209 type II AGNs at 0.3 < z < 0.83 from the 
Sloan Digital Sky Survey. The reliable stellar velocity dispersions (ct*) are obtained 
for 33 type II AGNs with significant stellar absorption features. According to the 
Lroiii] criterion of 3 x 10 8 L©, 20 of which can be classified as type II quasars. We 
use the formula of Greene & Ho to obtain the corrected stellar velocity dispersions 
(cr£). We also calculate the supermassive black holes masses from a% in these high- 
redshift type II AGNs. The [O III] luminosity is correlated with the black hole mass 
(although no correlation between the extinction-corrected [O III] luminosity and the 
black hole mass) , and no correlation is found between the Eddington ratio and the [O 
III] luminosity or the corrected [O III] luminosity. Three sets of two-component profile 
are used to fit multiple emission transitions ([O IIIJAA4959, 5007 and [O II]AA3727, 
3729) in these 33 stellar-light subtracted spectra. We also measure the gas velocity 
dispersion (a g ) from these multiple transitions, and find that the relation between 
<j g and <j% becomes much weaker at higher redshifts than in smaller redshifts. The 
distribution of < o g ja% > is 1.24 ± 0.76 for the core [O III] line and 1.20 ± 0.96 for 
the [O II] line, which suggests that a g of the core [O III] and [O II] lines can trace cr£ 
within about 0.1 dex in the logarithm of a%. For the secondary driver, we find that 
the deviation of <j g from er£ is correlated with the Eddington ratio. 
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1 INTRODUCTION 

Recent advances in the study of normal galaxies and active 
galactic nuclei (AGNs) are ample observational evidence for 
the existence of central supermassive black holes (SMBHs) 
and the relationship between SMBHs and bulge properties 
of host galaxies (Gebhardt et al. 2000; Ferrarese & Mer- 
ritt 2000; Tremaine et al. 2002; Begelman 2003; Shen et al. 
2005). We can use the stellar and/or gaseous dynamics to 
derive the SMBHs masses in nearby inactive galaxies. How- 
ever, it is much more difficult for the case of AGNs. With the 
broad emission lines from broad-line regions (BLRs)(e.g. H/3, 
Mg II, CIV; Ha), the reverberation mapping method and the 
empirical size-luminosity relation can be used to derive the 
virial SMBHs masses in AGNs (Kaspi et al. 2000; Vester- 
gaard 2002; McLure & Jarvis 2002; Wu et al. 2004; Greene 
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& Ho 2006a). It has been found that nearby galaxies and 
AGNs follow the same tight correlation between the central 
SMBHs masses (Mbh) and stellar bulge velocity dispersion 
(<r*) (the Mbh — <t* relation) (Nelson et al. 2001; Tremaine 
et al. 2002; Greene & Ho 2006a, 2006b), which also implied 
that the mass from reverberation mapping method is reli- 
able. 

According to AGNs unification model (e.g. Antonucci 
1993; Urry & Padovani 1995), AGNs can be classified into 
two classes depending on whether the central engine and 
BLRs are viewed directly (type I AGNs) or are obscured by 
circumnuclear medium (type II AGNs). In type I AGNs, by 
using the broad emission lines from BLRs (the reverberation 
mapping method or the empirical size- luminosity relation), 
we can derive virial SMBHs masses. It is not easy to study 
their host galaxies because their optical spectra are dom- 
inated by the non-stellar emission from the central AGNs 
activity. This is especially true for luminous AGNs, where 
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the continuum radiation from central source outshines the 
stellar light from the host galaxy. 

In type II AGNs, the obscuration of BLRs makes both 
the reverberation mapping method and the empirical size- 
luminosity relation infeasible to derive SMBHs masses. How- 
ever, we can use the well-known Mbh — o"* relation to de- 
rive SMBHs masses if we can accurately measure the stellar 
bulge velocity dispersion (cr*). There are mainly two dif- 
ferent techniques to measure cr*, one is the "Fourier- fitting" 
method (Sargent et al. 1977; Tonry & Davis 1979), the other 
is the "Direct-fitting" method (Rix & White 1992; Greene 
& Ho 2006b and reference therein). These years it has been 
successful to derive cr* through fitting the observed stellar 
absorption features, such as Ca H+K AA 3969, 3934, Mg 
IbAA 5167, 5173, 5184 triplet, and Ca IIAA 8498, 8542, 8662 
triplet, etc, with the combination of different stellar tem- 
plate spectra broadened by a Gaussian kernel (e.g. Kauff- 
mann et al. 2003; Cid Fernandes et al. 2004a; Greene & Ho 
2006b). 

On the other hand, Nelson & Whittle (1996) find that 
the gaseous velocity dispersion (cr g ) of [O III]A5007 from 
the narrow-line regions (NLRs) is nearly the same as cr, 
for a sample of 66 Seyfert galaxies, and suggest that the 
gaseous kinematics of NLRs be primarily governed by the 
bulge gravitational potential. Nelson (2001) find a relation 
between Mbh and o"[oiii] (the [O III]A5007 velocity disper- 
sion) for AGNs, very similar to the relation of Mbh — cr*, 
although with more scatter, which strongly suggests that a g 
can be used as a proxy for cr*. For lower-redshift type II 
AGNs with 0.02 < z < 0.3, Kauffmann et al. (2003) have 
investigated the properties of their hosts from the Sloan Dig- 
ital Sky Survey (SDSS) Data Release One (DR1), measured 
cr, and estimated the SMBHs masses from a* (Brinchmann 
et al. 2004). By using this sample, Greene & Ho (2005) have 
measured the gaseous velocity dispersion (cr g ) from multiple 
transitions ([O II] A3727, [O III] A5007, and [S II] AA6716, 
6731) and compared cr, and a g . They find that a g from 
these multiple transitions trace cr, very well, although some 
emission features they show considerable scatters. 

Type II quasars are the luminous analogs of low- 
luminosity type II AGNs (such as Seyfert 2 galaxies). The 
obscuration of BLRs makes quasars appear to be type II 
quasars (obscured quasars). Some methods have been used 
to discover type II quasars, but only a handful have been 
found. Recently, Zakamsa et al. (2003) present a sample of 
291 type II AGNs at redshifts 0.3 < Z < 0.83 from the SDSS 
spectroscopic data. About half are type II quasars if we use 
the [O III] A5007 line luminosity to represent the strength 
of the nuclear activity. What is the cr, — o g relation for type 
II quasars? And what are their SMBHs masses and the Ed- 
dington ratios Lboi/^Edd (i.e. the bolometric luminosity as 
a fraction of the Eddington luminosity)? 

Here we use the sample of Zakamsa et al. (2003) to study 
these questions for type II quasars. In section 2, we introduce 
the data and the analysis. Our results and discussion are 
given in Sec. 3. All of the cosmological calculations in this 
paper assume Ho = 70 km s _1 Mpc -1 , !Jm = 0.3, Qa = 0.7. 



2 DATA AND ANALYSIS 

With SDSS, Zakamsa et al. (2003) presented a sample of 291 
type II AGNs at redshifts 0.3 < z < 0.83. We downloaded 
these spectra from SDSS Data Release Four (DR4) and the 
spectra for 202 type II AGNs at redshifts 0.3 < z < 0.83 are 
obtained. SDSS spectra cover 3800-9200 A, with a resolution 
(A/AA) of 1800 < R < 2100 and sampling of 2.4 pixels per 
resolution element. The fibers in the SDSS spectroscopic sur- 
vey have a diameter of 3" on the sky, for our Type II AGNs 
sample at redshifts 0.3 < z < 0.83, the projected fiber aper- 
ture diameter typically contains about 90% of the total host 
galaxy light (Kauffmann & Heckman 2005), thus it is feasi- 
ble to observe significant stellar absorption features, which is 
the key point to accurately measure the stellar velocity dis- 
persion (cr*) in this paper. Since the SDSS spectra measure 
the light within a fixed aperture size, the estimated veloc- 
ity dispersions of more distant galaxies could be affected by 
the rotation of stars at larger physical radii than for nearby 
galaxies. In order to check the rotation contribution to line 
broadening, we examine the relation between the stellar ve- 
locity dispersion and the redshift, and find no correlation 
at all, which suggests that the line-broadening contribution 
from rotation of host galaxies is very small and negligible. 
The main reason is that cr* is derived through fitting heavy- 
element absorption lines (such as Mg IbA 5173, and Call K A 
3934, etc.), which are mainly from the old stellar population 
in the bulge. 

We first modelled the stellar contribution in the SDSS 
spectra of type II AGNs through the modified version of 
the stellar population synthesis code, STARLIGHT (ver- 
sion 2.0, Cid Fernandes et al. 2001; Cid Fernandes et al. 
2004a; Cid Fernandes et al. 2004b; Garcia-Rissman et al. 
2005), which adopted the new stellar library from Bruzual 
& Chariot (2003). The code does a search for the linear 
combination of Simple Stellar Populations (SSP) to match 
a given observed spectrum O x . The model spectrum M x is: 



M\(x, M\ , Ay, U*, cr*) = M\ 



.3=1 



®G(u*,a*)(l) 



where bj, A = L X SP (tj , Zj) / L x ^ p (tj, Zj) is the spectrum of 
the j th SSP normalized at A , r x = \Q-° A ^- A ^ is the 
reddening term, x is the population vector, M X() is the syn- 
thetic flux at the normalization wavelength, and G(«*, cr*) is 
the line-of-sight stellar velocity distribution, modelled as a 
Gaussian centered at velocity v* and broadened by cr*. The 
match between model and observed spectra is calculated by 
X 2 (x,M Xo ,A v ,v+,a i ,) = £x=iK°* - M x )w x ] 2 , where the 
weighted spectrum w x is defined as the inverse of the noise 
in O x . For more detail, please refer to Cid Fernandes et al. 
(2005). 

Prior to the synthesis, the Galactic extinction is com- 
puted by using the extinction law of Cardelli, Clayton & 
Mathis (1989) and the Ay value is taken from Schlegel, 
Finkbeiner & Davis (1998) as listed in the NASA/IPAC Ex- 
tragalactic Database (NED). The spectra are transformed 
into the rest frame defined by the redshift given in their 
SDSS FITS header. The spectrum is normalized at 4020A 
and the signal-to- noise ratio is measured in the S /N window 
between 4730 and 4780 A. Masks of 20-30 A around obvious 
emission lines are constructed for each object individually. 
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Because the redshift coverage of this type II AGNs sam- 
ple, we focus on the strongest stellar absorption features of 
Ca II K and the G-band, which are less affected by nearby 
emission lines. An additional /„ ~ power-law com- 

ponent (Watanabe et al. 2003) is used to account for the 
contribution from the scattered AGNs continuum emission, 
a traditional ingredient for modeling Seyfert galaxies since 
Koski (1978). Finally we check visually our spectral fitting 
results one by one. 

For our sample, the S/N in the S/N window varies be- 
tween 0.3 and 21.5. The fitting results for high S/N objects 
are usually better than those for low S/N ones. After in- 
specting the fitting results, we find that the fitting goodness 
(chi-square value) depends not only on the S/N (> 5 , in 
the given S/N spectral window), but also on the absorption 
lines equivalent widths (EW of Ca II K line > 1.5A). At last 
we select 33 type II AGNs, which had shown significant stel- 
lar absorption features and are well fitted to derive reliable 
measurements of stellar velocity dispersion. 

In order to check whether this sub-sample is represen- 
tative of the total sample of Zakamsa et al. (2003) with 
respect to the [O III]A5007 luminosity (£[0111] )> we plot the 
histograms of the Ljoiii] distribution for the sub-sample and 
total-sample (see Figure 1). Then T-test shows that these 
two populations are drawn from the same parent population 
with a possibility of 0.95. £[0111] is directly adopted from 
Table 1 in Zakamsa et al. (2003), the adopted raw I/[oiii] 
values and associated quantities are lower limits. Using the 
L[om] criterion of 3 x 10 8 L0(the logarithm is ~ 8.48), which 
corresponds to the intrinsic absolute magnitude Mb < — 23 
(Zakamsa et al. 2003), 20 objects can then be classified as 
type II quasars. Fig. 2 shows a fitting example of for SDSS 
J150117.96+545518.2 with S/N=20.5. The final results are 
presented in Table 1. All the fittings for 33 type II AGNs 
are appended in the Appendix. 

After subtracting the synthetic stellar components and 
the AGNs continuum, we obtain the clean pure emission- 
line spectra as shown in the top panel of Fig. 2, where we 
can analyze the pure emission-line profiles in detail by using 
the multi-component spectral fitting task SPECFIT (Kriss 
1994) in the IRAF-STSDAS package 1 . Because of the asym- 
metric profiles of the [O IIIJAA4959, 5007 lines, two sets of 
two-gaussian profiles are used in order to remove properly 
the asymmetric blue/red wings of the [O III] line. We take 
the same linewidth for each component, and fix the flux 
ratio of [O III]A4959 to [O III]A5007 to be 1:3. For three ob- 
jects, we can't fit the [O III] lines for their irregular [O III] 
lines (see Table 1). For the [O II] AA 3727, 3729 lines, we use 
two-gaussian profiles and fix their wavelength separation to 
the laboratory value; the ratio of the line intensities is al- 
lowed to vary during the fitting. The decomposition for [O 
II] lines is more difficult because of relatively low S/N and 
that the expected line widths are comparable to the pair 
separation (2.4A). For objects with available Ha, we also 
used one-gaussian profile to fit Nil AA 6548,6583 and Ha 
6563 lines. And Ha and Ha luminosities would be used to 
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do the extinction correction. For more details, please refer 
to Bian, Yuan & Zhao (2005, 2006). Our sample fitting for 
SDSS J150117.96+545518.2 is showed in Fig. 3. 



3 RESULTS AND DISCUSSION 

3.1 The uncertainties of the stellar velocity 
dispersion (cr,) and the gaseous velocity 
dispersion (<r 3 ) 

The derived stellar velocity dispersion was corrected by the 
instrumental resolutions of both the SDSS spectra and the 
STELIB library. Cid Fernandes et al. (2004a) have used their 
stellar population synthesis method to study a sample of 79 
nearby galaxies observable from the southern hemisphere, of 
which 65 are Seyfert 2 galaxies. The S/N in the S/N win- 
dow varies between 10 and 67. They compared their a* with 
values from the literature and found the agreement is good. 
They estimated that the uncertainty in <r„ is typically about 
20 km s _ . Recently, Cid Fernandes et al. (2005) applied 
their synthesis method to a larger sample of 50362 galax- 
ies from the SDSS Data Release 2 (DR2). Their derived a* 
is consistent very well with that of the MPA/JHU group 
(Kauffmann et al. 2003). The median of the difference be- 
tween the two estimates is just 9 km s _ . We have carefully 
checked our synthesis fitting result one by one and picked 
out 33 type II AGNs that are well fitted and the stellar ve- 
locity dispersion (<r*) are reliably derived. The spectral S/N 
for these objects are in the range of 5 to 21.5, most of which 
larger than 10, thus the typical uncertainty in a* should be 
around 20 km s~ . 

Recently, Greene & Ho (2006a, 2006b) used the direct- 
fitting method (Barth et al. 2002) to study the the system- 
atic biases of <r* from different regions around Ca II triplet, 
Mg lb triplet, and Call H+K stellar absorption features, 
which are introduced by both template mismatch and con- 
tamination from AGNs. They argue that the Ca II triplet 
provides the most reliable measurements of cr* and there 
is a systematic offset between ucaK and cr* derived from 
other spectral regions. For our higher-redhsift sample and 
the SDSS wavelength coverage 3800-9200A, it is impossi- 
ble to measure <r* from Ca II triplet. Therefore, for higher- 
redhsift type II AGNs, new observations around Ca II triplet 
are necessary in the future. Here we used the following for- 
mula to obtain the corrected velocity dispersion 0% (Greene 
& Ho 2006b), 

or* = (1.40±0.04)cr, -(71±5). (2) 

For three objects, <r* is near the instrumental resolution and 
the corrected crj is unreliable. These objects are excluded 
from further analysis, and denoted as f in Table 1. 

The gaseous velocity dispersion (cr g ) is obtained from 
full width half maximum (FWHM) of emission lines by 
assuming the Gaussian profile: crp OII ] = FWHM[OII]/2.35 
and af OIII] = FWHM[OIII]/2.35. Taking into account the 
SDSS spectrum resolution, the intrinsic a g derived from 
FWHMQOIII]) may be instrumentally broadened. The in- 
trinsic <T g value can be approximated by a g — (a 2 bs — 
(o"inst/(l + z)) 2 ) 1 / 2 , where z is the redshift. For the spectra 
from SDSS, the mean values of <7i n st are 74 km s _1 (the log- 
arithm is 1.87 dex) for [O II], and 60 km s~ 1 (the logarithm 
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is 1.78 dex) for [O III] (Greene & Ho 2005), respectively. 
The results of <7 g are listed in Table 1 (columns 7 and 8). 
After removing the effect of instrumental broadening, some 
objects become unresolved (see Table 1). Measurements of 
cr g below the resolution limit (74 and 60 km s~ 1 for [O II] and 
[O III], respectively) are not reliable. The error of er g is de- 
rived from the measurement error of the linewidth. For the 
linewidth, the typical error is about 10 per cent. However, 
the systematic effects are neglected, e.g. the uncertainties of 
the continuum subtraction and component decomposition 
(Bian, Yuan & Zhao 2005). 



3.2 The extinction correction of the [O III]A5007 
luminosity 

The [O III] line has the advantage of being strong and easy 
to detect in AGNs. The [O III]A5007 line is usually used as 
a tracer of AGNs activity (e.g. Heckman et al., 2004; Greene 
& Ho 2005). However, the [O IIIJA5007 luminosity is subject 
to extinction by interstellar dust in the host galaxy and in 
our Galaxy. Using SDSS data, Kauffmann et al (2003) found 
that for type II AGNs at z < 0.3 the mean value of the ratio 
of the corrected luminosity to the observational luminosity 

G°g( )i tne extinction factor at 5007A) is about 0.85 

dex and suggested that the extinction effect should be huge 
for type II quasars. The extinction correction is usually cor- 
rected by using the Balmer decrement, which is regarded as 
the best approximation. For 18 out of these 33 objects with 
available Ha measurements in the stellar-light subtracted 
spectra, we calculated the extinction correction of the [O 
IIIJA5007 luminosity, from the relation (Bassani et al. 1999) 



J [OIII] 



J [OIII][ 



r (H a /H 5 ) 



bs ,2.94 



(Ha/Hfl) 



(3) 



where an intrinsic Balmer decrement (H a /H^)o = 3.0 is 
adopted (Gu & Huang 2002). The range of the extinction 
factors at [OIII] is between -0.24 dex and 1.25 dex. The 
extinction-corrected [O III]A5007 luminosity is listed in Col. 
11 in Table 1. If we used the L[oiii] criterion of 3 x 10 8 L© 
(Zakamsa et al. 2003), 14 out of these 18 objects can then 
be classified as type II quasars. 

Using the SDSS AGN catalogue at MPA/JHU (Kauff- 
mann et al 2003) 2 , we obtained the extinction factor at 



5007A, log( 



). We tried to obtain a relation between 



J [OIIl] 
^[OIII 



and other observational parameters, such as er», 

L[oiii], Lh/3- However, no correlation is found. More careful 
work is required to deal with this problem in the future. 



3.3 The SMBHs masses and Eddington ratios in 
type II quasars 

Using the reverberation mapping method or the empiri- 
cal size-luminosity relation, it is impossible to estimate the 
SMBHs masses in type II quasars for the lack of emis- 
sion lines from BLRs. Here we use the formulae to de- 
rive the SMBHs masses in type II quasars from the stel- 
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lar velocity dispersion <r» (Tremaine et al. 2002), which is 
M BH (M ) = 10 8 - 13 K7(200 kms" 1 )) 4 02 . 

We calculate the Eddington ratio, Lboi/LEdd- We use 
the extinction-corrected [O III] luminosity (L? OIII ) as 
a surrogate for the AGNs luminosity (Heckman et al., 
2004; Greene & Ho 2005), L bo i = 3500Uf OIII] , to calcu- 
late the bolometric luminosity Lboi, and LEdd = 1-26 x 
10 38 Mbh/M© ergs s~ . Objects with no extinction cor- 
rected [O III] luminosity are discarded in the analysis rela- 
tive to the Eddington ratio. The results of SMBHs masses 
and Eddington ratios in type II AGNs are also presented in 
Table 1 (columns 13 and 14). We also calculate the SMBHs 
masses and Eddington ratios for the lower-redshifts type II 
AGNs at 0.02 < z < 0.3 presented by Kauffmann et al. 
(2003). For the typical uncertainties of 20km s _1 for a* = 200 
km s _1 , the errors of log cr* would be about 0.05 dex, cor- 
responding 0.17 dex for logMuH, and almost the same for 
Lboi/LEdd (Bian, Yuan, & Zhao 2005). Here we don't con- 
sider the cosmology evolution of Mbh— <t* relation (e.g. Woo 
et al. 2006), which is a question open to debate. 

In Fig. 4, we compare the distribution of SMBHs masses 
and Eddington ratios in the lower-redshifts sample and the 
higher-redshifts sample. It is found that the type II AGNs 
at higher-redshifts have higher SMBHs masses and higher 
Eddington ratios. 

It is well known that the Eddington ratio is an impor- 
tant parameter to describe the accretion process in AGNs. 
The [O III] luminosity is usually used as a surrogate for 
the AGNs luminosity (Heckman et al. 2004 and the ref- 
erence therein). In Fig. 5, we plot the [O III] luminosity 
and the corrected [O III] luminosity versus the SMBHs 
masses. Using a least-squares regression, we derive the cor- 
relation between Mbh and I/join] to be: log(L[ in]/LQ) = 
(6.83 ± 0.65) + (0.22 ± O.O8)log(M BH /M ). The correlation 
coefficient R is 0.45, with a probability of P„ u ii = 0.012 for 
rejecting the null hypothesis of no correlation. However, no 



correlation is found between Mbh and L 



[OIII] 



No correla- 



tion is found between the Eddington ratio and the [O III] 
luminosity or the corrected [O III] luminosity. From the pho- 
toionization model, the strength of [O III] is controlled by 
the NLRs covering factor, density, and ionization parame- 
ter (e.g. Baskin & Laro 2005). The relation between the [O 
III] luminosity and the AGNs bolometric luminosity is still 
a question to debate. 



3.4 The (T g — <7» relation 

The existence of a good correlation between stellar velocity 
dispersion (ct») and ionized gas velocity dispersion (cr g ) (e.g. 
Nelson & Whittle 1996) suggests that the gaseous kinemat- 
ics of NLRs in Seyfert galaxies be primarily dominated by 
the bulge gravitational potential, which is further confirmed 
by Nelson & Whittle (1996). Most recently, Greene & Ho 
(2005) have investigated a large and homogeneous sample 
of lower-redshift (0.02 < z < 0.3) Type II AGNs from the 
SDSS, and found that a s traces Though the gas kine- 
matics of NLRs are governed by the gravitational potential 
of the bulges, they also find that the accretion rate plays an 
important secondary role. 

Following Greene & Ho (2005), we study the a g — <r» 
relation for 33 Type II AGNs at redshifts 0.3 < z < 0.83 after 
measuring the gas velocity dispersion (cr g ) from the narrow 
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emission lines from NLRs, and the stellar velocity dispersion 
(at) from the Call H+K, G-band absorption feature, which 
is shown in Fig. 2. In Fig. 6, we showed the relation between 
a g and at- The correlation coefficients are 0.19, 0.27, 0.70, 
0.54 for figurs in Fig. 6, from left to right and from top 
to bottom. The relation between a g and at becomes much 
weaker at higher redshifts. We also found that this poor 
correlation is not due to the S/N. In order to qualify the 
comparison between a g and <rj, we calculate the distribution 
of < o-g/o-J >. The value is 1.24 ± 0.76 for the the core 
component of [O III] line, 1.20 ± 0.96 for the [O II] line. 
These suggest that a s of the the core component of [O III] 
and [O II] lines can trace at within 0.09 and 0.08 dex in the 
logarithm of at, respectively, and that the high-ionization 
[O III] line traces at as well as the low-ionization [O II] line. 
If we use the line width of [O III] core component to estimate 
the SMBHs masses from the Tremaine's Mbh — a relation, 
we would overestimate SMBHs masses by about 0.38 dex. 

For comparison, Greene & Ho (2005) derived the dis- 
tribution of < CT g /(T* > of lower-redshift type II AGNs with 
0.02 < z < 0.3, which are 1.34 ± 0.66 for the [O III] line, 
1.00±0.35 for the [O III] core line, and 1.13±0.38 for the [O 
II] line. We also calculate the the distribution of < a g /at > 
for the sample of Seyfert galaxies (Nelson & Whittle 1996); 
the value is 1.15 ±0.68, suggesting that a s of the [O III] line 
can trace at within 0.06 dex in the logarithm of at . Our 
results are thus consistent with theirs. 

To a first-order approximation, the line widths of the 
core component of [O III] and [O II] for both low-redshift 
and high-redshift Type II AGNs are primarily controlled by 
the gravitational potential of the bulges of host galaxies, 
and can approximately trace the stellar velocity dispersion. 
As we know, the errors in virial SMBHs masses derived from 
galaxy dynamics or size-luminosity relation is about 0.5 dex 
(e.g. Magorrian et al. 1998; Kaspi et al. 2000; Bian & Zhao 
2004a, 2004b). We also can use the line- width of the core 
component of [O III] or [O II] to estimate the black hole 
mass. 

In order to find the secondary effect of the line broad- 
ening in gas lines from NLRs for our higher-redshift narrow- 
line AGNs, we studied the relation between the deviation 
of crg from at (Act = loga g — logat) and the Eddington 
ratio (-Lboi/iEdd) (Greene & Ho 2005). Using the least- 
square regression, we find a median strong correlation be- 
tween these two dimensionless parameters. For the [O III] 
line, the relation is: Act = (-0.035 ± 0.057) + (0.202 ± 
0.058)£og(L bo i/L Edd ), R=0.71, P nuU = 0.00441. For the [O 
II] line, the relation is Act = (0.047 ± 0.058) + (0.179 ± 
0.059)Zo g (L bol /L Edd ), R=0.65, P nuU = 0.0091 (see Fig. 7). 
We also find the comparable correlation between Act and 
Mbh- These results confirm that the nuclei accretion pro- 
cess and/or nuclei SMBHs would affect the linewidth of gas 
lines from NLRs, although the primary driver is the gravi- 
tational potential of the bulge. 



4 CONCLUSION 

The stellar population synthesis code is used to model the 
stellar contribution for a sample of 209 type II AGNs at 
redshifts 0.3 < z < 0.83 from SDSS. According to the L[oiii] 



criterion of 3 x10 s L , 20 can be classified as type II quasars. 
The main conclusions can be summarized as follows. 

• The reliable a, are measured for 33 type II AGNs with 
significant stellar absorption features. We use the formula 
of Greene & Ho to obtain the corrected stellar velocity dis- 
persions (at and SMBHs masses are calculated from the 
Mbh — o% relation. A median strong relation between the 
[O III] luminosity and the SMBH mass is found (although 
no correlation between the extinction-corrected [O III] lu- 
minosity and the SMBH mass); no correlation is found be- 
tween the Eddington ratio and the [O III] luminosity or the 
extinction-corrected [O III] luminosity. 

• The gas velocity dispersion (<r 9 ) in NLRs is measured 
using three sets of two-gaussian profiles to fit [O III]AA4959, 
5007 and [O II]AA3727, 3729) in these 33 stellar-light sub- 
tracted spectra. We find that the relation between a g and 
at becomes much weaker at higher redshifts. 

• The distribution of < a g /at > is 1.24 ± 0.76 for the 
core [O III] line and 1.20 ± 0.96 for the [O II] line, which 
suggests that a g can trace at within about 0.1 dex in the 
logarithm of at. The deviation of a g from at is correlated 
with the Eddington ratio. 
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Table 1. Results for type II quasars, sorted by the value of S/N. The columns are as follows: (1) name; (2) redshift; (3) uncorrected 
stellar velocity dispersion in units of km s" 1 ; (4) FWHM of [O III] A 5007 A in units of km s _1 ; (5) FWHM of [O II] A 3727 A in units 
of km s ; (6) log of corrected stellar velocity dispersion by equation 2 (in units of km s _1 ); (7) log of the [O III] a in units of km s —1 ; 
(8) log of the [O II] a in units of km s — , (9) log of [O III] luminosity in unit of Lq; (10) log of [O II] luminosity in unit of Lq; (11) log 
of extinction-corrected [O II] luminosity in unit of Lq; (12) signal-to-noise ratio (S/N); (13) log of the black hole mass in units of solar 
mass; (14) log of the Eddington ratio, j: cr* is near the instrumental resolution and the corrected u% is unreliable, which is excluded from 
future analysis. 



8 W. Bian, Q. Gu, Y. Zhao, L. Chao, Q. Cui 



50 



40 - 



30 - 



20 - 



10 - 







-i — | — i — | — i — | — i — | — i — | — i — | — i — | — i — | — i — | — i — r 



-i 1 r- 



J i I i I lL 



sub-sample 
total sample 



' ■ I i L 




i ■ i 



5.5 6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 

'°9 ( L [ o„, ] /L sun) 



Figure 1. The distribution of £[0111] f° r total sample and the sub-sample. 
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Figure 2. Sample fit of the synthetic population model for SDSS J150117. 96+545518. 2. Top: the observed spectra (top black curve, 
shifted up for clarity), the synthetic spectra (middle red curve), and the residual spectrum (bottom black curve). Bottom: the region 
around Ca H+K AA 3969, 3934 and G-band (left); the region around Mg IbAA 5167, 5173, 5184 triplet (right). 
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Figure 3. Sample multi-component fitting of the [O II] AA3726,3729 and [O III]AA4959, 5007 lines for SDSS J150117.96+545518.2: 
modeled composite profile (thick solid red line), individual components (the dotted green lines), the residual spectrum (lower panel). 
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Figure 5. The [O III] luminosity (left) and the corrected [O III] luminosity (right) plotted against the SMBHs masses derived from the 
stellar velocity dispersion a. Solid line is the best fit. 



12 W. Bian, Q. Gu, Y. Zhao, L. Chao, Q. Cm 





log o. (km s" ) 



log cy, {km s ) 






1 + 1 v- 







log J" ""(km s" 1 ) 



Figure 6. Top left panel is <r g versus <rj for the [O III] line; top right panel is <r g versus crj for the [O II] line; bottom left is cr g 0111 ^ 
versus a g OII ^\ bottom right is <r g versus cr* for the the sample of Nelson & Whittle (1996). The uncertainty in a is typically 20 km s _1 . 
The solid line denotes 1:1. 




Figure 7. Act = logcr g — logcr% plotted against the Eddington ratio, Lboll T^Edd- Left panel is for the [O III] line and right panel is for 
the [O II] line. Solid lines show the best fits. 



